On-site Coulomb interaction and the magnetism of (GaMn)N and (GaMn)As. 
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We use the local density approximation (LDA) and LDA+U schemes to study the magnetism of 
(GaMn) As and (GaMn)N for a number of Mn concentrations and varying number of holes. We show 
that for both systems and both calculational schemes the presence of holes is crucial for establishing 
ferromagnetism. For both systems, the introduction of U increases derealization of the holes and, 
simultaneously, decreases the p-d interaction. Since these two trends exert opposite influences on 
the Mn-Mn exchange interaction the character of the variation of the Curie temperature (Tc) 
cannot be predicted without direct calculation. We show that the variation of Tc is different for 
two systems. For low Mn concentrations we obtain the tendency to increasing Tc in the case of 
(GaMn)N whereas an opposite tendency to decreasing Tc is obtained for (GaMn)As. We reveal 
the origin of this difference by inspecting the properties of the densities of states and holes for 
both systems. The main body of calculations is performed within a supercell approach. The Curie 
temperatures calculated within the coherent potential approximation to atomic disorder are reported 
for comparison. Both approaches give similar qualitative behavior. The results of calculations are 
related to the experimental data. 
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I. INTRODUCTION 

After the discovery of the ferromagnetism of 
Gao.947Mno.o53As with the Curie temperature as high 
as llOKi the diluted magnetic semiconductors (DMS) 
became the subject of intensive experimental and theo- 
retical studies being considered as promising materials 
for semiconductor spintronics. Although much new un- 
derstanding of the physics of DMS has been won by the 
investigation of (GaMn)As, the Curie temperature of this 
prototype system could not be raised above 150Ki£ The 
design of technologically useful DMS materials with the 
Curie temperature exceeding the room temperature re- 
mains a challenge. 

The theoretical prediction by Dietl et al^ of the high- 
Tc ferromagnetism in (GaMn)N played an important 
role in the formulation of the directions of the search 
for new ferromagnetic DMSs. A number of experimental 
studies indeed detected the ferromagnetism of (GaMn)N 
samples with the Curie temperature higher than the 
room temperature^^ However, the experimental data 
on the magnetic state of (GaMn)N are strongly scattered 
and range from a paramagnetic ground state to the fer- 
romagnetic state with a very high Curie temperature of 
940 K. 

New proof of the complexity of the magnetism of 
(GaMn)N was given by the recent magnetic circular 
dichroism measurements by Andoi The measurements 
were performed on a high-Tp sample and led Ando to 
the conclusion that the (GaMn)N phase in this sample 
is paramagnetic. The ferromagnetism of the sample was 
claimed to come from an unidentified phase. On the other 
hand, in a recent preprint by Giraud et a\m a high temper- 
ature ferromagnetism was detected in the samples with 
a low Mn concentration of about 2%. The authors rule 
out the presence of precipitates in the system and argue 
that the room temperature ferromagnetism is an intrin- 
sic property of (GaMn)N. Thus the situation is contro- 



versial and calls for further research efforts. The purpose 
of the present work is a theoretical comparative study of 
(GaMn)As and (GaMn)N with a focus on the influence 
of the on-site Coulomb interaction on the Curie tempera- 
ture of both systems. The calculations are performed on 
the basis of the density functional theory (DFT) within 
both LDA and LDA+US approach. 

Numerous LDA calculations of the electronic struc- 
ture of (GaMn) As have been reported (see, e.g., reviews 
Refsi^Siii and more recent calculations RefsiiSii&i^ii^) . 
For (GaMn)N the number of such calculations is still 
small (see RefsiiSiiLi2iiS and recent preprint 20 ). The 
study of the properties of both systems which goes be- 
yond the LDA is at the very beginning. There are a 
short report on the LDA+U calculation for (GaMn)As 21 
and a very recent article^ on LDA+U calculation for 
both systems. We mention also the preprint on the self- 
interaction corrected LDA (SIC-LDA) calculation^ for 
these systems. 

The present work is novel in a number of aspects. First, 
we investigate the characteristics of the holes and discuss 
them from the viewpoint of the mediation of the ferro- 
magnetic exchange interaction between the Mn impuri- 
ties. Second, we calculate interatomic exchange inter- 
actions and the Curie temperatures. Third, we separate 
the contributions of the antifcrromagnetic superexchangc 
and the ferromagnetic kinetic exchange to the Curie tem- 
peratures. In all studies we focus on the influence of the 
on-site Coulomb interaction. 



II. CALCULATIONAL SCHEME 

The calculational scheme is discussed in RefsiiS^ to 
which the reader is referred for more details. The scheme 
is based on the DFT calculations for supercells of semi- 
conductor crystals with one Ga atom replaced by a Mn 
atom. The size of the supercell determines the Mn con- 



2 



centration. To calculate the interatomic exchange inter- 
actions we use the frozen-magnon technique and map the 
results of the total energy of the helical magnetic config- 
urations onto a classical Hciscnberg Hamiltonian 



H, 



eff 



5> 



(1) 



where Jij is an exchange interaction between two Mn 
sites and ej is the unit vector pointing in the direc- 
tion of the magnetic moment at site i. 

To estimate the parameters of the Mn-Mn exchange 
interaction we perform calculation for the frozen-magnon 
configurations: 



const, 



q R t 



(2) 



where Qi and 4>i are the polar and azimuthal angles of 
vector e^, Ri is the position of the ith Mn atom. The di- 
rections of the induced moments in the atomic spheres of 
the atoms of the semiconductor matrix are kept parallel 
to the z axis. 

It can be shown that within the Heisenberg model 
the energy of such configurations can be represented in 
the form 



E(6,q) = E (9) -sin 2 6»J(q) 



(3) 



where Eq does not depend on q and J(q) is the Fourier 
transform of the parameters of the exchange interaction 
between pairs of Mn atoms: 



J (q) = J( y ex p(*q ■ R oj) 
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(4) 



Performing back Fourier transformation we obtain 
the parameters of the exchange interaction between Mn 
atoms: 

q 

The Curie temperature is estimated in the mean-field 
(MF) approximation 
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We use a rigid band approach to calculate the exchange 
parameters and Curie temperature for different electron 
occupations. We assume that the electron structure cal- 
culated for a DMS with a given concentration of the 3d 
impurity is basically preserved in the presence of defects. 
The main difference is in the occupation of the bands 
and, respectively, in the position of the Fermi level. 

In the main body of the LDA+U calculations we 
use U — 0.3Ry that corresponds to the value deter- 
mined experimentally 2 ^, and is used in some previous 
calculationsi 19 i 21 i 24 The dependence of the Curie tem- 
perature on U is illustrated at the end of the paper. The 
results of the supercell calculations are compared with 
the results of the calculations within the coherent poten- 
tial approximation (CPA). 2 - 5 
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FIG. 1: (Color online.) The spin-resolved total DOS of 
(GaMn)As and (GaMn)N for Mn concentration of 6.25%. 
Calculations are performed within LDA and LDA+U ap- 
proaches. For (GaMn)N, both zinc-blende (ZB) and wurzite 
(WU) structures are presented. The spin-up/down DOS is 
shown above/below the abscissa axis. The total DOS is given 
per chemical unit cell of the semiconductor. 



III. CALCULATIONAL RESULTS AND 
DISCUSSION 

A. Density of states 

We begin the discussion with the properties of the den- 
sity of states (DOS). In Figs. UHwe present the DOS for 
x = 6.25%. The main features of the DOS discussed be- 
low are common for all Mn concentrations studied. The 
replacement of one Ga atom by a Mn atom adds to the 
energy structure five valence bands related to the Mn3d 
spin-up states. On the other hand, there are only four 
valence electrons added. As a result, both (GaMn)As 
and (GaMn)N have one hole per Mn atom. (GaMn)N 
is half-metallic in both LDA and LDA+U calculations 
with only majority-spin states present at the Fermi level. 
(GaMn)As is half-metallic in LDA calculation and pos- 
sesses small contribution of the minority-spin states in 
LDA+U calculation. In all cases the systems are char- 
acterized by high spin-polarization of the electron states 
at the Fermi energy. This property is in good agreement 
with a recent measurement of high carrier spin polariza- 
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FIG. 2: (Color online.) The same as in Fig. 0but for partial 
Mn3d-DOS. 



tion in (GaMn)AsiS 6 . 

Despite this similarity, important differences between 
(GaMn)As and (GaMn)N are obtained in both LDA and 
LDA+U calculations. In LDA, the spin-up impurity 
band of (GaMn)As merges with the valence band. On 
the other hand, in (GaMn)N the impurity band lies in 
the semiconducting gap. This is valid for both wurzite 
and zinc-blende crystal structures of (GaMn)N. These 
features are in agreement with the results of previous 
calculations, 10 ! 11 ! 12 ! 16 ! 17 ! 18 ! 19 ! 20 In all LDA calculations 
the impurity bands have large Mn3d contribution (Fig. 
0). 

In the LDA+U calculations, the impurity band of 
(GaMn)As disappears from the energy region at the top 
of the valence band. Almost all spin-up Mn3d states 
lie at about — 0.4Ry below the Fermi level. In contrast, 
(GaMn)N still possesses impurity band which lies now 
closer to the valence band. The Mn3d contribution to 
the impurity band decreases compared with the LDA- 
DOS but is still large (Fig. [TJ). 



B. Properties of the holes 

Now we turn to the discussion of the properties of the 
holes. We again present the results for one Mn concen- 
tration (Fig. |3J) and begin the discussion with the spa- 
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FIG. 3: (Color online.) Distribution of the hole in LDA and 
LDA+U calculations. Numbers in parentheses give the num- 
ber of atoms in the coordination sphere. The hole part is 
given for one atom. The distance from Mn is given in units 
of the lattice parameter. The narrow rectangulars to the left 
and right of the ordinate axis present the Mn3d contribution 
into the hole for LDA and LDA+U calculations. 



tial distribution of the hole in the LDA calculations. In 
(GaMn)As 29% of the hole is located on the Mn atom. 
The first and second coordination spheres of As include 
four atoms each and contain 35% and 14% of the hole re- 
spectively. The rest is distributed between Ga and empty 
spheres. 

In (GaMn)N the holes are much stronger localized 
about the Mn atom. For the ZB structure we find 50% 
of the hole in the Mn sphere. First and second coordi- 
nation spheres of N contain respectively 33% and 5% of 
the hole. Although the part of the hole located at the 
distant atoms decreases compared with (GaMn)As it is 
still sizable. This property allows the mediation of the 
exchange interaction between the Mn atoms. Similar re- 
sult we obtain for the wurzite (GaMn)N (not shown). 

In the LDA+U calculations, the holes become more 
delocalized (Fig. [3Jl. In (GaMn)As, the Mn contribution 
decreases to 9%. On the other hand, the contribution 
on the first and second coordination spheres of As in- 
creased to 45% and 23%. The contribution of the Ga 
atoms changes much less. For (GaMn)N in ZB struc- 
ture the corresponding numbers are 25%, 48% and 12%. 
Again similar behavior is obtained for (GaMn)N in the 
wurzite structure. 

The increased derealization of the holes is a factor 
favorable for the mediation of the exchange interaction 
between the Mn atoms and therefore for a higher Curie 
temperature. There is however another factor that is 
equally important for the establishing the long-range fer- 
romagnetic order. This factor is the strength of the p-d 
exchange interaction which is the physical origin of the 
interaction mediated by the holes. This second factor 
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is less easy to characterize quantitatively. A convenient 
quantity for such a characterization is the p-d exchange 
parameter usually denoted as N(3, In the mean-field ap- 
proximation N (3 describes the spin-splitting of the semi- 
conductor valence-band states which appears as a conse- 
quence of the interaction of these states with Mn3d elec- 
trons. However, our calculations (Fig. EJ) show that the 
spatial form of the hole states can differ strongly from the 
unperturbed semiconductor states making the mean-field 
definition of the p-d interaction poorly founded. 

The experimental estimations of N/3 made for 
(GaMn)As on the basis of different experimental tech- 
niques vary strongly from large values of \N/3\ = 3.3 cV 27 
and N/3 = 2.5 eV— to much smaller values N/3 = —1.2 
e\& and \N/3\ = 0.6 eV. 29 (We use the traditional sign 
convention: negative N/3 corresponds to antiparallel di- 
rections of the spins of the d and p states.) A possible 
reason for this variation is strong p-d hybridization in 
III-V DMS. Note that for many II- VI systems N/3 is well 
defined experimentally 30 and allows good DFT descrip- 
tion within the LDA+U approach— In the present paper 
we do not consider the problem of the estimation of the 
N/3 parameter in systems with strong p-d hybridization. 
For characterization of the strength of the p-d interac- 
tion we use the value of the Mn3d contribution to the 
hole (Fig. EJ. 

For both systems the LDA+U calculations give a 
strong decrease of the Mn3d contribution to the hole 
(Fig. EJ. For (GaMn)As, the Mn3d contribution drops 
from 16% to a small value of 4%. For (GaMn)N in ZB 
structure we get 43% and 19% and for (GaMn)N in WU 
structure 41% and 16% (the results for WU structure are 
not shown in Fig. EJ). Important that in (GaMn)N the 
LDA+U the hole still contains large Mn3d contribution. 

The presence of two competing trends in the properties 
of holes - increasing derealization, on the one hand, and 
decreasing p-d interaction, on the other hand - makes 
direct calculation of Tc necessary to predict the influence 
of the Hubbard U on the Curie temperature. 



C. Curie temperature 

In Fig. El we show the Curie temperature for a number 
of concentrations of Mn impurities— 

The calculated Curie temperatures (Fig. EJ show again 
strong difference between (GaMn)As and (GaMn)NiS In 
(GaMn)N, the Curie temperature increases with account 
for on-site Coulomb interaction for all Mn concentrations 
considered. On the other hand, in (GaMn)As the value 
of Tc decreases for concentration 3.125%, 6.25% and 
12.5%. 

For comparison we present (Fig. EJ the Curie temper- 
ature calculated within the coherent potential approxi- 
mation. The supercell approach employed in the present 
paper and the CPA provide two opposite limits in the 
treatment of the atomic disorder: The supercell calcula- 
tion deals with an ordered atomic pattern. On the other 
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FIG. 4: (Color.) Tc as a function of the Mn concentration. 
The superexchange contribution Ts and hole contribution 
Tj? are shown (see the text for definitions of these contri- 
butions). The calculations are performed for ZB structure. 
For (GaMn)N, the results for WU structure for x = 12.5% 
are presented for comparison. The lower part of the figure 
shows the Curie temperature calculated within CPA. L+U is 
abbreviation for LDA+U. 



hand, the CPA assumes complete atomic disorder with 
the effects of the short range order neglected. Important 
that the influence of Hubbard U on the Curie tempera- 
ture obtained within CPA is qualitatively similar to that 
obtained within the supercell calculations: In (GaMn)As, 
U gives a lower Tc for small x and a higher Tc for larger 
x. In (GaMn)N, account for Hubbard U leads to an in- 
crease of Tc for the whole Mn concentration range stud- 
ied. 

To get a deeper insight into the formation of the Curie 
temperature we study the dependence of Tc on the band 
occupation. In Fig. El we present the results for two Mn 
concentrations of 3.125% and 12.5%. Here, n = corre- 
sponds to the nominal number of carriers (one hole per 
Mn atom). The kink at n = 1 marks the point where the 
valence band is full and the conduction band is empty. 
In Fig. El we show the dependence of the leading inter- 
atomic exchange parameter— as a function of the band 
occupation. Indeed, the behavior of the Curie tempera- 
ture (Fig. EJ repeats in gross features the behavior of the 
exchange parameter (Fig. EJ. In particular, the decrease 
of the number of holes below the nominal value of one 
hole per Mn atom leads to the decrease and change of the 
sign of the exchange parameter. Thus our calculations 
show the crucial role of the valence-band holes for medi- 
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FIG. 5: (Color online.) The Curie temperature as a func- 
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for, respectively, £=3.125% and x=12.5% as a function of the 
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ating ferromagnetism in both (GaMn)As and (GaMn)N. 
This result is obtained in both LDA and LDA+U ap- 
proaches and agrees with the commonly accepted pic- 
ture of the ferromagnetism in DMS used in the model- 
Hamiltonian studies (see, e.g., review o 33l34 n and references 
therein) . 

It is useful to present the Curie temperature in the 
form Tc=Ts+T# . Here Tg is the Curie temperature at 
n = 1 and specifies the exchange interaction mediated by 
the completely filled valence band. This contribution can 
be related to the Anderson's superexchange ^ T# is the 
difference between the Curie temperatures at n = 1 and 
n = and specifies the contribution of the hole states. 

The value of T5 is negative for all concentrations in 
both LDA and LDA+U approaches reflecting the prop- 
erty that the exchange interaction through completely 
filled bands is antiferromagnetic. The hole contribution 
Th is always positive and is responsible for the ferromag- 
netism of the system. 

The analysis of Fig. 0] shows that the contribution 
of the antiferromagnetic superexchange is in general not 
small and must be taken into account in the estimation 
of the Curie temperature. The absolute value of the anti- 
ferromagnetic superexchange is always larger in the LDA 
case. The decrease of the superexchange for nonzero U 
can be explained by an increased energy distance between 
the occupied and empty Mn3d states. 

The influence of Hubbard U on the hole contribution 
Tjf differs strongly for two systems. In (GaMn)N, LDA 
and LDA+U give very similar results up to x=12.5%. For 
larger x the LDA value of Tc is higher. This shows that 
for low Mn concentrations the increased derealization 
of the holes is compensated by the decreased p-d inter- 
action. For larger x the increase of the p-d interaction 
becomes more important since the distance between the 
Mn atoms decreases. For (GaMn)As, the strong decrease 
of the p-d interaction under the influence of Hubbard U 
is not compensated by the increased hole derealization 
and the value of Tjj drops strongly for small x. 

Comparison of the results of the calculations with 
experiment shows that in the case of (GaMn)As the 
account for Hubbard U does not make the agreement 
with experiment better. Indeed, the calculated Curie 
temperature becomes too low. A possible reason for 
this can be an underestimation of the Mn3d contribu- 
tion into the hole states within the LDA+U calculation: 
these states disappear from the Fermi-level region al- 
most completely (Fig.[Q. The dynamical mean-field the- 
ory (DMFT) which treats the on-site correlations beyond 
the LDA+U can improve the agreement with experiment 
since the DMFT provides an account for the quasiparti- 
cle states^ that are neglected within the LDA+U. Before 
the LDA+DMFT study is performed the decision about 
more appropriate method for the study of the ferromag- 
netism in (GaMn)As should be made in favor of the LDA 
scheme which gives good correlation with experiment^ 

In (GaMn)N, the Mn3d states are present at the Fermi 
level because of the strong hybridization with the N2p 
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FIG. 7: (Color online.) The Curie temperature of (GaMn)As 
with 3.125% of Mn as a function of the band occupation for a 
number of U values. In the insert, the Tc value as a function 
of U for nominal electron number n = 0. 

states. We obtained an increase of Tc which is in cor- 
relation with the prediction by Dietl et al£ and with 
some of the experimental results. Our LDA+U results 
for (GaMn)N correlate in many respects with the picture 
proposed recently by Dietl et al2I in the discussion of the 
Zhang- Rice limit for this system. They conclude that the 
high Tc can appear in (GaMn)N as a result of strong p-d 
interaction within the Zhang-Rice polaron which is effi- 
ciently transfered due to the relatively large extent of the 
p-wave functions^ 

Figure [7| illustrates the dependence of the Curie tem- 
perature on the value of U for (GaMn)As with 3.125% 
of Mn. We obtain a systematic trend of decreasing scale 
of the variation of Tc with variation of the band occu- 
pation. The insert shows Tc(U) for nominal electron 
number n = 0. The dependence is monotonous though 
non-linear. 

Within our LDA+U calculation we did not obtain the 
orbital ordering found in the SIC calculation^ We, how- 
ever, do not exclude that such a state can be stabilized 
within the LDA+U approach. Our numerical experiment 
with the impurity band split "by hand" to simulate the 
effect reported in Rcf. 20 indeed resulted in a nonferro- 
magnetic ground state. 

IV. CONCLUSIONS 

We use the local density approximation (LDA) and 
LDA+U schemes to study the magnetism of (GaMn)As 
and (GaMn)N. We show that for both systems and both 
calculational schemes the presence of holes is crucial for 



establishing ferromagnetism. The introduction of U in- 
creases derealization of the holes and, simultaneously, 
decreases the p-d interaction. Since these two trends ex- 
ert opposite influences on the Mn-Mn exchange interac- 
tion the character of the variation of the Curie tempera- 
ture (Tc) cannot be predicted without direct calculation. 
We show that the variation of Tc is very different for two 
systems. For low Mn concentrations we obtain the ten- 
dency to increasing Tc in the case of (GaMn)N whereas 
an opposite tendency to decreasing Tc is obtained for 
(GaMn)As. We reveal the origin of this difference by 
inspecting the properties of the densities of states and 
holes for both systems. The main body of calculations is 
performed within a supercell approach. The Curie tem- 
peratures calculated within the coherent potential ap- 
proximation to atomic disorder are reported for compar- 
ison. Both approaches give similar qualitative behavior. 
We make a contact between calculational results and the 
strong scattering of the experimental data on Tc of the 
systems. 

Note that the present experimental and theoretical sit- 
uation does not allow to give a unique answer on the ques- 
tion which of two schemes, LDA or LDA+U, provides 
better description of the diluted magnetic semiconduc- 
tors. This answer can be different for different systems 
and even samples. Our previous calculations for the II- VI 
DMS (ZnMn)Se2£ gave strong arguments for the superi- 
ority of the LDA+U approach in this case. In the present 
study of (GaMn)As, the LDA+U method with a moder- 
ate U underestimates the strength of the p-d interaction 
and leads to too low Curie temperature for low Mn con- 
centrations. For (GaMn)N, the experimental situation is 
too unclear to allow for a well-founded decision concern- 
ing the superiority of the theoretical scheme. The com- 
parative study reported here aims to contribute to both 
the increase of the understanding of the systems studied 
and the formulation of the theoretical approaches most 
suitable for the investigation of concrete DMSs. The ef- 
forts must, however, be continued and we hope that the 
results of this paper will be stimulating for further stud- 
ies. 
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